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ABSTRACT 

We make a detailed analysis of the BBp{DDp) and B*Bp{D*Dp) strong couplings gBBpidDDp) 
and gB*Bp{gD*Dp) using QCD light cone sum rules(LCSR). The existing some negligence is 
pointed out in the previous LCSR calculation on gB*Bp{9D*Dp) and an updated estimate is 
presented. Our findings can be used to understand the behavior of the B, D ^ p semileptonic 
form factors at large momentum transitions. 
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1. INTRODUCTION 

At the present time, there is an increasing interest in exclusive B decays in order to explore 
the sources of CP violation. However, the definite interpretations for the relevant experimental 
data demand that we have the ability to precisely compute the physical amplitudes. It is 
nonperturbative QCD dynamics not being dealt rigorously with that would hinder us from 
doing such a desired calculation. One is forced to use some approximate methods. Lattice 
QCD simulation is the most trustworthy approach to nonperturbative QCD effects, with no 
parameters or assumptions, but the precision of calculation is limited by the available computing 
resources and some certain restriction exists[|I|, for example, in describing b ^ uor s transitions. 
The formulation of QCD factorization formula^ is based on the first principle and is viewed 
as a great progress in phenomenology of heavy fiavors; however, the underlying long distance 
effects included in a series of the hadronic matrix elements still confront us. Heavy quark 
symmetry can very well apply to describing the systems including one heavy quark, but is of 
less predictive power for heavy-to-light transitions. Although QCD sum rule method]^ builds 
its underlying physical assumptions on the field theory, and exhibits its decided superiority 
in dealing with some of nonperturbative quantities, such as decay constants, hadronic matrix 
elements and strong coupling constants, a problem with it is that the resulting form factors for 
heavy-to-light transitions can not behave very well in the heavy quark limit mq oo. The 
successes of perturbative QCD in treating numerous hard exclusive processes have excited the 
occurrence of many excellent works. The most prominent of them is the development of QCD 
light-cone sum rules (LCSR)[|, It is the striking advantage of this approach to describe 
heavy-to-light transitions in a way consistent with the universally accepted physical picture 
that nonperturbative QCD dynamics occupies an dominant place and perturbative hard gluon 
exchanges contribute only a subleading effect in that case. LCSR approach takes the basic 
correlator, in which the proper current operators are sandwiched between the vacuum and an 
on shell light meson state, and adopts the operator product expansion (OPE) around the light 
cone ~ instead of at the small distance x ~ 0. Nonlocal matrix elements occurring in 
this approach, which encode all the information on large distance dynamics, are parametrized 
in terms of a set of so-called light-cone wavefunctions classified by twist, which describe the 
momentum distributions of the quarks inside the relevant light mesons. This is identical with 
a summation over all the condensate terms in the short distance OPE. Consequently, the 
resulting form factors for heavy-to-light transitions exhibit the correct behavior with heavy 
quark mass, averting the problem with the traditional sum rules. This method has extensively 
been accepted and has found its successful applications, such as the investigations of the form 
factors for heavy-to-light transitions at small and intermediate momentum transfers [5- 11] and 
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of the strong couplings between heavy and hght mesonsP, 12, 13|, since its presentation in 



Very recently, it has been generalized to study the nonfactorizable effects in B ^ TrvrfQ and 



to probe heavy-to-light form factors in the whole kinematically accessible ranges [p!5| . The 
technical details of LCSR can be found, for instance, in Q], while for a detailed comparison 
with traditional sum rules, see 0. 

An investigation of strong interactions between heavy mesons and a light vector meson is 
of important phenomeno logical interest, and especially the relevant strong couplings can be 
employed to understand the behavior of the corresponding heavy-to-light form factors at large 
momentum transfer in a pole dominance model, which has proven to be exact in an effective chi- 
ral Lagrangian approach [P!^. The off-shell B*Bp{D*Dp) strong coupling gB*Bp{gD*Dp), in fact, 
has been calculated in the standard LCSR approach |TB[; however, it calls for a reexamination 
due to some negligence existing in calculations. Another nonperturbative quantity deserving 
of investigation is the BBp{DDp) strong coupling gBBpidDDp), which turns out to be equally 
important. Motivated by all these facts, in this paper we make a systematic study on them. 

This presentation is organized as follows. The following Section is devoted to a detailed 
derivation of the sum rules for gBBpigoDp) and gB*Bp{gD*Dp), using LCSR method. Then we 
give a numerical analysis of the resulting sum rules, including a discussion of error estimates, 
in Sec. 3. The last Section give to a simple summary. 

2. LCSR'S FOR THE STRONG COUPLINGS 

The BBp and B*Bp strong couplings gBBp and gB*Bp can be defined as 

{p {q, e) B (p) \B {p + q)) = gBBpe\x) ■ P, (1) 

(p (g, e) B*{p, r^)\B {p + q)) = -gB*Bpepa(3,V*'q''efx)P^ . (2) 

The resulting findings can easily converted into the corresponding c-quark meson cases. Accord- 
ing to the general strategy of QCD sum rules, it is needed to construct an adequate correlator in 
order to obtain an qualitative estimate for gBBp and gB*Bp- As usual, we use two correlators of 
the following forms as the starting points of LCSR calculations on gBBp and gB*Bp, respectively, 



F (p, (p + q)) = ij d^xe'P"" (p (g, e) \Tu (x) i-f^b (x) , b (0) ij^d (0) 







F{p',{p + q)')e*-p, (3) 



Gpip,p + q) = rf'xe*^-"(p(g,e)|Tll(x)7^6(a;),6(0)^75C^(0)|0; 

= G{pMp + qf)e,^f,,q-e*^p\ (4) 
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For the correlator (3), isolating the pole contribution of the lowest 0~ B meson and parametriz- 
ing these from the higher 0~ states in a form of dispersion integral, the hadronic form of the 
invariant function F ^p^, (p + q)^^ may be written as 



{rriB-p^) {rriB-p + q) (si-p^) S2-{p + q) 

with a double integral starting from the same threshold parameter sq, which should be set in 
the neighborhood of the squared mass of the first excited 0" B meson. Similarly we have for 
the invariant function {p, {p + q)), 

qH r 2 ^ + )2\ ^ rnB*m%fBfB*9B*Bp + f f {si, S2) dsids2 



nib {m%* - p2) [777,^ - (p + qy] J J {si - p^) [s2 - (p + q)"^] 

QCD calculations of the underlying correlators may be allowed, on the other side, for the 
negative and large values of p^ and (p + ?)^, in which case the b quarks travel only a small 
distance x and therefore the operator product expansion (OPE) goes effectively in powers 
of the deviation from the hght cone x'^ ^ 0. We would like to work in the case where the 
interactions of the b quarks with the background field gluons are omitted, since their influence 
on the sum rules is, as always, negligibly small. On contracting the b quark operators into a 
free propagator, 

(0|r6WM0)|)^pl^/dUe-l±^, (7) 

one gets, 

FQ^-^{p',{p + qf) = (2^// ciWA;^^^^e'(^-'=)-[^'^(p(g,e)|rK(a;)7^ci(0)|0) 

-mb{p(q,e)\Tu{x)d{0)\0)]. (8) 

The nonlocal matrix elements {p {q, e) \Tu {x) 7^0? (0) | 0) and {p{q, e) \Tu {x) d (0) | 0) define the 
light cone wavefunctions of the p meson[17-19] as, 

q ■ X 

+ (eW* - £ due^-^-g^^ («, /.) 

ip-xf ' Jo '^'''^ I ' 

(p(g, e) \Tu (x) d (0)| 0) = -i/2fjml{e^ ■ x) £ diie-^-/i^ (li, ptb). (10) 

/p stands for the usual vector decay constant of the p meson and fj is defined as ^0|Jo'pj,'u|p^ = 
i/J (e^i^^qu — e^^q^y, <^|| {u, /J,) is the leading twist-2 wavefunction, g^^^ {u, n) and /i|(M, /x^) refer 



{p{q,e)\u{x)j,d{0)\0) = f^mA^^qJ due^^^'^ 

I q ' X t/ 
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to the twist-3 ones, and both A {u, fi) and C {u, /i) have twist-4, which parametrize the mass 



corrections. (f\\ {u,fi), {u,fi) and h^^{u,fib) are normahzed as Jq duf (u) = 1, while C {u,fi) 
satisfies Jq duC (u) = 0. A tedious but straightforward calculation yields 



du- 



.u 



ml — {p + uq) m\ 
h\f\u,iJb) 



(p + q) _ 



+fjmlmb j du— ^ 

° mi — {p + uq) 



-- fo'm^n [ du I ■ 

r'^ p Jo y 



ml — {p + uqY 



+ 



ml 



ml — {p + uqY 



X 



A {u) + 8C (u) 



fill 



In deriving Eq.(ll), it proves to be convenient for a partial integration to introduce the auxiliary 
functions / (u) = Jq f (v) dv for f{u) = ip\\{u, jj), g^^\u, jj) and A{u, fj), and C{u) = C{v)dv, 
C (u) = Jq C (f ) dv, which is equal to zero when u = 1. The twist-3 contribution from g^j^\u, fi) 
vanishes exactly due to cancellations in the partial integrations. 

Furthermore, it is indispensable to convert Eq.(ll) into a form of dispersion integral for 
upcoming continuum substraction. The relevant QCD spectral density can easily be obtained 



by virtue of the technique suggested in [^]. In the following we consider only the twist-2 and 
-3 terms and give a detailed deviation. First of all, we perform a double Borel transformation 



= —p^ — ^ M^, Q\ = — {p + qf M|, for the twist-2 term ( the term proportional to the 



inverse ml — {p + qf disappears after doing that ), obtaining 



b{mIqI)b{mIqI) [ du 



M2 



with Mo = Mf / (Mf + M|) and M 



ml — {p + uq) 

2 M?M2 



7^K+™'"«(^-"o^](^ll (uo), (12) 



. The symmetry of the correlator makes it natural 



to set Ml = M| so that the wavefunctions (u) may take its value at the symmetric point 
Uq = 1/2. Further, making a replacement Mf l/cri,M| l/o"2 in Eq.(12) yields 



B 



0-2 



Q 



du- 



ml — {p + uqY 



f (cri,0-2) • 



(13) 



Finally, we take the function f {01,02) 
mation in the variables ai and a2'- cri 
density reads 



((Ji, (J2) and perform one more Borel transfor- 
1/si and a"2 l/s2- The resulting QCD spectral 



(QCD) 
Ptw2 



Sl,S2j 



S1S2 
fp^p'PW 



0-1 ) 5 — ,(T2) / (o-i,a2 



S2 



5 {Sl - S2) ( S 



1 2 



6 S2 - m, 



1 2 



.(14) 



Applying all the same procedure to the twist-3 part, we have 

pSS^\suS2) = fjrnlrnb(ml + ^ml^ (^^^ 6 (^si - ml - ^mfj 6 (^S2 - ml - ^mfj 

(15) 

which means that the twist-3 part receives no continuum substraction. With Eqs. (14) and 
(15), F^^^ [p^^ {p + q)^^ can be expressed as 



6 {si - S2) (si - m,l - \mfj Q (s2 - ml - ^ 



4^1 



fjmlmhhl (i) 



Si -p2) U2 - {p + qY 



(ml + \ml - p^) (ml + \ml - (p + qf) 



du 



+ 



mt 



f^l ~ (P + '*^?) fTT'b ~ (P + 



A{u) + 8C{u) . 

(16) 



Making the Borel improvement p"^ Mf, {p + q)"^ M| for both the theoretical and 
hadronic expressions, which suppresses the higher state and twist-4 contributions, and then 
equating them by the use of the quark-hadron duality ansatz, the final sum rule for f^gBBp 
reads. 



2 2 

flgBBp = —e'^ \ fpmpM' 
m^ 



_ "0 

e 



X 



^mbfJh['^{-]-f,m, 



.,1).8C(1 



1 + 



m^ 
M2 



(17) 



Now let's turn from this topic to a discussion of the sum rule for the B*Bp strong couphng 
gB*Bp by expanding the relevant correlator (4) around the light cone — 0. Utilizing Eq.(7) 
it follows immediately that 

1 



G^^^^ipMp + qf) 



I I d'xd'k , e'(P-^)-"[r(p|rTZ(x) 7^7.75^(0)10) 



(2vr) 

+m, (p|T«(a;) 7^75^(0)10)] 



(18) 



The 7 algebraic relation 7,^71/75 = —^ei^^a/s'^'^^ +91x^15 can help us to write it down in a preferred 
form. 



G<^^^^pMp + qf) 



{2n) 
1 



d xd hj 



ml — k"^ 



i{p—k)x 



[k^ {p\u{x)-i^d{Q)\Q) 



e^uapk" (p u {x) a'^'^d (0) o) + m^ {p \Tu (x) 7^75^ (0) | 0) 



.(19) 
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The nonlocal matrix element (p \u{x) 75c/ (0)| 0) is exactly vanishing, as required by the parity 
conservation of strong interactions. The light cone expansions of the other two matrix elements 
reads [ffn, O] respectively, 



P 



uix) a 



'"d (0) 







*/9 a 







^(A) ■ ^ 2 



due''"'-'' 
1 



V5± {u) + —rripX At {u) 



(q-x) 







and 



X'' 



1 



*/3 c 



q ■ X Jo 



due''"'-''CT {u) 

Jo 



{p \Tu (x) 7^75^ (0)1 0) = ^fpTUpe^af^^q^'e^x"' 



due''"'''g^^\u,fib). 



(20) 



(21) 



(p± (u) and g^_^^ (u) stand for the twist-2 and -3 wavef unctions, respectively, and obey f (u) du = 
1; the others are associated with twist-4 operators and parametrize the p mass corrections, 
among which both Bt (u) and Ct (u) abide by f (u) du = as it stands. 

As with the gsBp case, in practical calculations we use the definitions Bt (u) = Jq Bt (f ) dv, 
Bt {u) = £Bt (v) dv and Gt (u) = /g" Gt (v) dv. Substituting Eqs. (20) and (21) into (19) we 
gain the theoretical expression for G^'-'^{p'^, {p + q)'^), 



G^''^{pMp+qf) = fj{l 



du- 



ip± (u) 



{p + uq) 



2 + 



duATiu) 



2ml 



ml — {p + uq) 
1 





ml — {p + uqY 



+ -mhmpfp 



du- 



9V{.u) 



(22) 



ml — {p + uqf' 

It is interesting to note that the p mass effect is only offered by the twist-4 wavefunction At{u). 
Obviously, the QCD spectral densities concerning the leading and next-to-leading terms are, 
except for a constant factor, the same as those in the Qbbp case, respectively. Omitting details, 
we end up with the following sum rule for the product fB*fBgB*Bp, 



fB*fBgB*Bp 



rrib 



m^B-TTlB 



+l^p 



2 I 2 
2a7^ 



ml 



e 



4™pJ — e 



MV± ( 2 



M2 



e 



\+4"'p)At - 



I f -^(ml + ^ml) (a) /I 

+^fpmbmpe t 4 p)gy ( _ 



(23) 



It should be understood that the Borel variables Mf and M| have been taken equal once again, 
for the B and B* mesons are nearly degenerate in mass. This enables the relevant wavefunctions 
to take values at m = 1/2, to high accuracy. 



We would like to emphasize that the previous LCSR results for gB*Bp and (yf^tDplTSl are 
questionable, because of a missing factor 1/2 in front of the term —e^aXp'^'^^ in writing down 
the 7 algebraic relation •y^'jalb = —^/'^^iJ.aXpO'^'^+giJ.a, and a misuse of the subtraction procedure 
e~M^(™'!'"^4™'p) e~ AJ^l^b+i^v) _ e~i^ imposed on the twist-3 parts. 

3. Numerical Results 

The parameters to need fixing for a numerical estimate are those concerning the B* or 
B mesons and those describing the p meson. The former contain the decay constants fs 
and Jb*, mass parameters tub, mB* and rub, and threshold parameter sq. We use rub = 
4.8 ±0.1 GeV, TUB = 5.279 GeV and niB* = 5.325 GeV. With the two point sum rules 
formulated in0, the values of the decay constants fs and fB* are fixed at fs = 115 MeV 
and fs* = 125 MeV, corresponding to nib = 4.8 GeV, sq = 33 GeV^ and the leading order 
in as- Also, all the parameters of the p meson involved become now numerically available. 
We take as inputs the experimental values nip = 770 MeV, fp = 198 ± 7 MeV and the QCD 
sum rule result /J = 152 ± 9 MeV flSl at the scale Ub = 



m 



B 



2.5 GeV. Concerning 

the light cone wavefunctions appearing in our sum rules, there have been many discussions on 
them in the literature. It is in ||l^ that QCD sum rule method is first applied to study the 
twist-2 distribution amplitudes of vector mesons. Later on, a more systematic discussion was 
||18|| , where results of [1^ were critically examined and updated. Very recently, the 



given m 



authors of ||T9[ took further the meson mass corrections into consideration by introducing some 
higher twist distributions in the light cone expansions of the relevant nonlocal matrix elements, 
extending the work [|18[ by an additional use of the QCD equations of motion. The yielded 
findings, some of which will be used in our numerical analysis, have found applications in 
phenomenology of exclusive semileptonic and radiative B decays. The explicit forms of the 
light cone wavefunctions in relation to our sum rule calculations are 



h'^'^ (m, pb) 
(m, Pb) 
At {u, Pb) 
A {u,pb) 

C{u,pb) 

with the coefficient (pb) = 



6u (1 -u)(^l + 4 {pb) ^ [5 (2n - 1)' 

Qu (1 -u){l + 4 {Ph) \ [5 (2u - 1)' - 
6n(l - n)(l + 0.15[5(2n - 1)^ - 1]), 
6u(1-m)(1 + [5(2m-1)=^-1]), 



(24) 



4 20 8 ■ 



30^2(1 -u)^ 



3 10 10 ■ 

2 + y + y 

0.17 ±0.09 and a| {pb) 
8 



uHl-u)\ 



0.16 ± 0.09, ^3 = 0.023 and ^4 = 0.13. 



Having all the input parameters at hand, we could carry out the numerical calculations. It 
is a critical step towards deriving a reliable sum rule prediction to look for a reasonable range 
of the Borel parameters. The standard procedure requires that the terms proportional to the 
highest inverse power of the Borel parameters stay reasonably small, which can fix the lower 
limit of the fiducial Borel interval, and that the higher resonance and continuum contribution 
should not become too large, which may determine the upper limit of the allowed range. For 
the two sum rules in consideration, we find that the Borel intervals to satisfy the above criteria 
are respectively 8 < < 14 GeV^ for the gsBp case and 8 < < 15 GeV^ for the gB*Bp 
case, where the twist-4 wavefunctions contribute less than 6% and 5% and the high states at 
the orders lower than 23% and 25%, respectively. The figure 1 shows the sensitivity of the sum 
rules for f^gsBp and fB*fBgB*Bp to the Borel parameters. Prom the corresponding sum rule 
"windows", we arrive at flgsBp = 0.071±0.002 GeV^ and fB*fBgB*Bp = 0.082±0.004 GeV, the 
uncertainties quoted being due to the variations of M^. By means of values of decay constants 
obtained previously it is immediate to get the desired sum rules for the strong couplings gBBp 
and gB*Bp- If using the central values for all the relevant sum rule results, we have gBBp — 5.37 
and gB*Bp = 5.70 GeV~^. 

For a better understanding of the overall uncertainties in the coupling constants, it is highly 
advisable to employ the analytic forms instead of the numerical results in Eqs.(17) and (23) 
for the decay constants. When rrib keeps fixed while Sq changes between 32 — 34 GeV"^, the 
resulting variations relative to the central values amount to ±6% for gBBp and to ±8% for 
gB*Bp, if the uncertainties due to the Borel parameters are included. The influences on the sum 
rules can be investigated of the uncertainty in m^, by considering a correlated variation of 
and So in the individually allowed ranges. Requiring the strong couplings gBBp and gB*Bp to 
take values only if the sum rules for the relevant decay constants show the best stability, we 
observe that the induced changes are typically of orders 7% and 6% respectively. It is of course 
important to investigate further the uncertainties from the wavefunctions. The simplest way to 
test the sensitivity of the sum rules to model wavefunctions is by putting all the corresponding 
nonasymptotic coefficients to zero, namely by using their asymptotic forms which are mode- 
independent and completely dictated by perturbative QCD. The resulting sum rules deviate 
from their individual central values by about 7% in the gBBp case and by about 8% in the gB*Bp 
case. The effects observed in such a way come from an extreme treatment and therefore are 
anyway being overestimated. Taking it into account that the twist-4 wavefunctions bring only 
a correction of about 4% to both sum rules for gBBp and gB*Bp, we can reasonably conjecture 
that the uncertainties due to neglected yet higher twists would be at best of the same orders as 
the twist-4 corrections. At present, the total uncertainties in the sum rules for gBBp and gB*Bp 
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can conservatively be estimated to be about 25% and 27%, respectively, by adding linearly up 
all the considered errors. 

The same procedures may be used for a numerical discussion of LCSR for goDp and gD*Dp- 
The relevant parameters are taken as rric — 1.3 GeV, itid — 1-87 GeV, mo* — 2.01 GeV, 
fo = 170 MeV, fo* = 240 MeV and sq = 6 GeV^ In addition, we have to evolute the 
wavefunctions to a lower scale /Xc = — m"^- Using the standard criteria the fiducial 

intervals of turns out to be 4 < < 8 GeV^ in the goDp case and 5 < < 8 GeV^ 
in the go* Dp case. The stability of the sum rules for f})gDDp and fD*fDgD*Dp is illustrated in 
Fig. 2. Wc have fj^goDp = 0.11 GcV^ and fD*fDgD*Dp = 0.17 GeV, with the negligibly small 
uncertainties due to M^. As for the strong couplings goDp and gD*Dp, the resulting sum rules 
are predicted to be goDp — 3.81 and go* Dp — 4.17 GeV~^, the total uncertainties being about 
23% and 25%, respectively. 

4. SUMMARY 

We have made an intensive study on QCD interactions between heavy mesons and a light 
vector meson within the framework of LCSR. A detailed deviation of the sum rules is presented 
for the relevant strong coupling constants gBBpiguDp) and gB*Bp{gD*Dp) and a systematic nu- 
merical analysis, including a painstaking investigation of the uncertainty arising from all the 
possible sources of error, is made. An existing negligence is pointed out in the previous LCSR 
calculation on gB*Bp and go* Dp, and an updated LCSR result is formulated. 

The obtained predictions can be used to estimate the couphngs for the other charge states 
using the relations from isospin symmetry. Also, it is straightforward to investigate the BgBK*, 
B*BK* and B*BsK* strong couplings, and the corresponding those in c quark meson case by 
making a corresponding parameter replacement in the relevant sum rules formulated. 

The numerical results presented here should be updated, once our understanding of the 
meson wavefunctions, 6-quark mass and decay constants became more clear, and the QCD 
radiative corrections are included in sum rule calculation. 
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FIGURE CAPTIONS 



Fig.l: The stability of LCSR for the products f^gBBpiFig.l (a)) and fB*fBgB*Bp{Pig.l (b)), 
with rUb = 4.8 GeV and Sq = 33 GeV^. 

Fig. 2: The stabihty of LCSR for the products /D25'DDp(Fig-2(a)) and /D*/D5'£)*Dp(Fig-2(b)), 
with nic — l.S GeV and sq — 6 GeV^. 
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